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Abstract 

Burgers equation is a quasilinear partial differential equation, proposed in 1930's to model 
the evolution of turbulent fluid motion, which can be linearized to the heat equation via the 
celebrated Cole-Hopf transformation. This paper introduces and studies in detail general ver- 
sions of stochastic Burgers equation with random coefficients, in both forward and backward 
sense. Concerning the former, the Cole-Hopf transformation still applies and we reduce a for- 
ward stochastic Burgers equation to a forward stochastic heat equation that can be treated 
in a "pathwise" manner. In case of deterministic coefficients, we obtain a probabilistic repre- 
sentation of the Cole-Hopf transformation by associating the backward Burgers equation with 
a system of forward-backward stochastic differential equations. Returning to random coeffi- 
cients, we exploit this representation in order to establish a stochastic version of the Cole-Hopf 
transformation. This generalized transformation allows us to find solutions to a backward 
stochastic Burgers equation through a backward stochastic heat equation, subject to addi- 
tional constraints that reflect the presence of randomness in the coefficients. In both settings, 
forward and backward, stochastic Feynman-Kac formulae are derived for the solutions of the re- 
spective stochastic Burgers equations, as well. Finally, applications that illustrate the obtained 
results are presented to a controllability problem and to a pricing/hedging problem arising from 
mathematical finance. 
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1 Introduction 



Burgers equation plays a very important role in the theory of differential equations and applied 
mathematics. It is a quasilinear parabolic partial differential equation (PDE) of the form 



where U(t, x) is an unknown velocity field, which is to be determined by the initial condition 
U(0,x) = uo(x), and v is a viscosity term; hereafter explicit dependence of the subsequent fields 
is suppressed for ease of notation and subscripts will denote partial derivatives with respect to 
the corresponding variable. This equation can be considered either in the whole real line (x € R) 
or in a bounded interval with Dirichlet or Neumann boundary conditions. It has been proposed 
in the 1930's by the Dutch scientist J.M. Burgers as a simple model for the dynamics of the 
Navier-Stokes equations in one spatial dimension [U [6]. As such it has been used to model the 
dynamics of one dimensional pressureless turbulence in fluid flows and the hope was that it would 
help to understand a lot of the intricate structure of this fascinating subject. One of the great 
breakthroughs in this study was the discovery by J. D. Cole [10] and E. Hopf [26] in the 1950's of 
a transformation that reduces this equation to the heat equation, thus allowing the derivation of 
exact solutions in closed form. In particular by the transformation U = — In V to a new variable 
V(t,x), where we take without loss of generality that v = 1, the Burgers equation transforms to 
the linear heat equation 



This transformation, which has been named the Cole-Hopf transformation,!], allows us to construct 
explicit analytic solutions for the Burgers equation. 

This breakthrough ended temporarily the career of Burgers equation as a modeling tool for tur- 
bulence, however, it turned it into a scientific paradigm (a) of a nonlinear PDE that is linearizable 
with the use of a simple tranformation (in some sense the PDE analogue of the Ricatti equation) 
and (b) of a benchmark model that can be used to understand the basic features of the interaction 
between nonlinearity and dissipation. The scientific community soon found many new and inter- 
esting applications for the Burgers equation, other than its initial ones in fluid mechanics. For 
instance, as a model for condensed matter physics, in statistical physics (a continuous time version 
of ballistic deposition models within the framework of the Kardar-Parisi-Zhang (KPZ) model), 
in cosmology within what is referred to as the Zel'dovich approximation [43j, in the modelling of 
traffic flow [9] , and even in economic theory [24\ [25j H2] . 

However, the Burgers equation developed a parallel, theoretical, and abstract mathematical 
existence beyond its dominant presence in applications. Motivated by the intention to reinstate 
Burgers equation as a model for turbulence, the nonlinear dynamics community turned its attention 
to the randomly forced Burgers equation. This introduced a class of quasilinear stochastic PDEs 

It is worth noting that a version of Burgers equation was proposed in 1906 by Forsyth, who in fact proposed a 
version of the linearising transformation, and a similar equation was further proposed by Bateman in 1915. For the 
history of the Burgers equation and the Cole-Hopf transformation see [19] . 



u t + uu x = - 
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(SPDEs) of the form 



U t + UU x = -U xx + F(t,x,u) 



where now F(-,x, •) is a stochastic process that acts as the forcing term. This turns the solution 
into a random field, whose properties are to be determined by the data of the problem, i.e., the 
forcing term and the initial condition. Depending on the properties of the random forcing term, 
this equation has to be treated accordingly, i.e., pathwise or in the ltd sense 



in the case where the forcing term F(t,x,oj) can be modelled as an infinite dimensional Brownian 
motion. 

The introduction of randomness in Burgers equation produced a number of very interesting 
new directions; directions connected with dynamical systems aspects of the equation, e.g. existence 
and properties of invariant measures (see for instance the important contributions by E et al [IT] 
or Goldys and Maslowski |21j). directions related to various questions on the well-posedeness of 
the equation in various functional settings using techniques from infinite dimensional stochastic 
analysis (see for instance the important contributions of Da Prato, Deboussche, Nualart and others 
|14[ \15\ [23"1 [22j [30l [37] ) , interesting connections with geometry (see for instance the contributions 
by Cruzeiro and Malliavin |12j or Davies, Truman, and Zhao [16J), connections with the theory of 
superprocesses in [3] etc. This theoretical work was inspired by issues related to turbulence (see 
e.g. [U [39l HO]) but also led to many new exciting applications (e.g. the work of Kiefer [29J that 
connects the random Burgers equation with polymer models, see also the review paper of Bee and 
Khanin [2 J on Burgers turbulence and references therein for details on other possible applications). 

Another interesting problem related to the Burgers equation, both for the deterministic as well 
as for the stochastic version, is that of its optimal control or controllability. This is an important 
problem both from the point of view of theoretical considerations as from the point of view of 
applications. Temam and coworkers [8] studied this problem from the viewpoint of application in 
turbulence control and gave some positive answers to the question of whether a feedback control 
law can be a feasible way to drive the fluid flow to a desired final state. Such questions were 
further elaborated upon by Da Prato and Debussche [13] with the use of dynamic programming 
techniques. However, there has been a renewed interest on the issue of controllability for Burgers 
equation; see for instance [TJ ETJ [20] for recent related work. These problems bring up the following 
two interesting directions. 

The first one is related to the well-posedeness of the forward stochastic Burgers equation, 
complemented with the feedback control law predicted by the dynamic programming equation. 
Following |13] such a closed loop equation may take the form 



where $(£,£/) is a functional specifying the feedback control law. The exact form of <5 depends 
on the cost criterion which is to be optimized, and at this point is irrelevant to our discussion. 





(3) 
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What is important, though, is that it introduces the need to study more complicated forms of the 
forward stochastic Burgers equation, that includes a potential like term $>(t,U). This term may 
well be a random field itself. 

The second one is related to the issue of controllability. When one considers the problem of 
whether it is feasible to find a control procedure W(t, U) which drives the system to a desired final 
state £ at a given time T, then it is appropriate to look at the original system as a final value 
problem with U(T,x) = £, rather than as an initial value problem. For instance, if one wishes to 
find W so as to drive the system 

U t -UU x = ~U xx + W(t,U) 

to the desired state U (T, x) = £ then one should study the well-posedeness of the resulting backward 
equation in the sense of the existence of a random field U and a mapping W(t, U) that satisfy 
the above problem. There are close connections between the problem of controllability and the 
well-posedeness of this backward problem; see for instance the discussion of the famous Hilbert 
Uniqueness Method (HUM) proposed by J. L. Lions |32j and J. M. Coron for the control of 
linear and nonlinear spatially extended systems, respectively. This powerful idea requires a major 
revision so that it may be applicable to SPDEs, on account of problems related to the adaptivity of 
the solution to the filtration generated by the noise process when time is reversed... Such questions 
led to the development of the theory of backward SPDEs (BSPDEs) (see for instance the work by 
Ma and Yong in [35J and [36 J that have recently found applications in |38| , [1] or |18| to "pathwise" 
utility optimization problems emerging from mathematical finance. Using this powerful theory we 
may consider the problem of controllability of the stochastic Burgers equation, by addressing the 
problem of well-posedeness of a backward stochastic Burgers equation that has the form 

dU= (-^U XX -UU X + <S>(U,Z)\ dt + ZdW(t), 
U(T,x)=£, 

where £ is a given random final condition, and now we need to look for the pair of random fields 
(U,Z). As expected the unknown random field Z is related to the control procedure needed to 
drive the system to the desired final state. 

In all the above, the stochasticity was supposed to be an additive term, imposed in the system 
through the external driving force F(t,x,u). However, the randomness may well be inserted in 
the model in terms of random coefficients. For instance, one may consider a model which is of 
the form of the Burgers equation, but with the viscosity term v being a random field rather than 
a constant parameter. Such a situation may be a good model for fluid turbulence, in which the 
viscosity term is determined through the distribution of turbulent eddies in the flow. Another 
case where a Burgers equation with random coefficients could arise, is when we consider a random 
potential. This may either arise as a physical model (e.g. as a model for a conducting fluid in 
the presence of external random electromagnetic fields) or in terms of a controlled problem, in 
the spirit of (J3j) - Other possible motivations could arise in the context of economic theory, and in 
particular within the context of portfolio selection in market models with random coefficients. 

It is therefore the aim of this paper to study forward as well as backward stochastic Burg- 
ers equations with random coefficients. Such problems have not been studied extensively in the 
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literature. One of our main concerns is to study whether a generalization of the Cole-Hopf trans- 
formation can be obtained, that allows us to transform the quasilinear Burgers equation with 
random coefficients to a linear heat equation. If this holds, then we may find solutions via this 
approach for the Burgers equation possibly in closed form, a fact that will allow us to obtain 
interesting information on their properties. Such information will provide important input for the 
relevant models. On the other hand, even if the solution of the heat equation in analytic form is 
no more possible, due to the presence of the random coefficients, the reduction to a linear equation 
will allow us to deduce interesting qualitative information on the solutions of the Burger equation, 
such as comparison principles, Feynman-Kac representations etc. 

In the forward case, we employ the Cole-Hopf transformation to a forward version of stochastic 
Burgers equation with random coefficients, linearizing it to a stochastic heat equation that does 
not contain a stochastic integral term; that is, it can be solved pathwise as a deterministic one. 
A special version of forward stochastic Burgers equation with constant coefficients was studied 
in |41j . The backward case is far more delicate and our method of approach is to associate the 
stochastic Burgers equation with a system of ordinary forward-backward stochastic differential 
equations (FBSDEs) with random coefficients, through the four step scheme introduced in [33J; 
an exposition of the theory on FBSDEs may also be found in [35]. To do so, firstly we establish 
this relation both for the Burgers equation and for its corresponding heat equation in the case 
of deterministic coefficients, developing a probabilistic approach to the Cole-Hopf transformation 
that allows us to reinterpret it as a point transformation between the two associated FBSDE 
systems. Using then this approach, we generalize this transformation to backward stochastic 
Burgers equations with random coefficients. This generalization is by no means trivial and reflects 
interesting features concerning the nature of the equation. Through this extended transformation, 
we find the most general version of backward stochastic Burgers equation that is linearizable and 
reducible to a stochastic heat equation, subject to additional constraint equations that now appear 
due to the randomness of the coefficients. Examples of completely solvable backward stochastic 
Burgers equations are also presented for particular families of random coefficients. 

In both cases, forward or backward, the resulting stochastic heat equation is linear and much 
easier to handle. Thus from the linear system and the generalized Cole-Hopf transformation we 
construct solutions to stochastic Burgers equations with random coefficients and obtain stochastic 
Feynman-Kac type representations for them. Our results indicate that the backward stochastic 
Burgers equations, whose study has been neglected so far in the literature, have more intricate 
structure than their forward counterparts and can lead to interesting applications. In particular, 
the controllability of a backward stochastic Burgers system that drives it to a predetermined final 
state amounts to its solvability subject to a suitably selected initial state control. On the other 
hand, in a tax regulated financial market with a money market and a stock, a small investor may 
price and hedge a contingent claim, whose future value depends on the volatility of the stock, by 
selecting his portfolio according to the solution of an FBSDE system and, in turn, to the solution 
of the associated backward stochastic Burgers equation. 

A summary of the paper is as follows. In Section [2] we introduce a forward and a backward 
general version of stochastic Burgers equation with random coefficients, and present the necessary 
mathematical preliminaries. Section [3] connects the backward deterministic Burgers equation with 
a class of FBSDEs, via the four step scheme, developing a probabilistic approach to the celebrated 
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Cole-Hopf transformation. In Sections 2] and [5] we cope with the linearization of the already intro- 
duced forward and backward stochastic Burgers equations with random coefficients, respectively. 
Stochastic Feynman-Kac formulae are established for their solutions in Section [6l Finally, Section 
[7] illustrates the obtained results with applications to controllability and mathematical finance. 

2 Random Burgers Equations: Introduction of Two Models and 
Mathematical Preliminaries 

The main goal of this section is two-fold. The first goal is to introduce two versions of stochastic 
Burgers equations with random coefficients, a forward and a backward one, and motivate them as 
generalizations of the deterministic Burgers equation in the presence of noise and randomness in 
the coefficients. The second goal is to set the necessary notation and functional framework, and 
also present briefly the necessary mathematical tools that will be used in the paper. 



2.1 Two versions of the random Burgers equation 

In contrast with the deterministic case where the forward and the backward problem in time can be 
related with the use of the simple time inversion transformation t — > — tp, in the stochastic case with 
or without random coefficients this is no longer true, on account of technical but vital for the nature 
of the problem difficulties related with the adaptedness of the solution to the filtration generated 
by the noise process. Therefore, in the stochastic case the forward and the backward problems have 
to be formulated and treated differently. In both problems, we consider a 1-dimensional Brownian 
motion W(-) on some filtered probability space (Q,J-,P;¥^ with F = J-(-) being the natural 
filtration generated by W(-), augmented by all the P-null sets in T. Generalizations to higher 
dimensional or even infinite dimensional noise are feasible but are not pursued in the present work. 
In addition, a, a, g, b, e, s, m, and / are assumed to be square integrable, real-valued random 
fields defined on [0, T] x R x for a terminal time T > 0, such that for fixed x S R they are 
F-progressively measurable. 

We consider first the forward version of the stochastic Burgers equation with initial condition 
given by a square integrable random field p : R x J7 — > R; that is, 



dU 



- a 2 (t,x)U xx + a(t,x)UU x +g(t,x)U^>(t,U)+b(t,x)U x 



+ e(t, x)U + s(t, x)-$ x (t, U) + m(t, x)#(t, U) + f(t, x) 

+ ^(t,U)dW(t), 0<t<T, xGR, 
U(0,x)=p(x), xGR. 



dt (4) 



Here ^f(t, U) is a known function (or even random field) of U, pre-described by the form of the 
model. It reflects the effect of the fluctuations on the evolution of U, i.e., how the effect of W(-) 



2 eventhough the parabolic nature of the equation may introduce certain difficulties related to the unique contin- 
uation of the solution for negative times 
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transfers to an effect on U. In the simple case where ^(t, U) is a constant, we obtain a stochastic 
Burgers equation with additive noise. 

The motivation for equation (JH) is rather straightforward. One may consider any physical 
situation where the evolution of the state of the system in question is modeled in terms of Burgers 
equation. Then assume that this system is subject to external sources, whose behavior is random 
and subject to fluctuations (this accounts for the Wiener process term in the above system), while 
at the same time the coefficients of the model are also subject to uncertainties. This uncertainty is 
introduced by the random coefficients, in a sense that the coefficients of the model are functions, 
rather than constants, of a random process that evolves in parallel with the state of the system. 
This random process depends on the external noise term, therefore we assume that the random 
coefficients are functions which are adapted to the filtration generated by the Wiener process. 

To illustrate the above points let us consider the example of turbulence modeling: the Wiener 
process models a body force for the fluid which presents fluctuations around an average body 
force. The random coefficients would correspond to a randomly changing viscosity term in the 
fluid. This is a very reasonable modeling assumption since for example the effective viscosity of 
the fluid depends on the eddie's formed in the turbulent flow which in turn depend on the random 
fluid velocity, which in turns depends on the external random body force. The above discussion 
can be easily transferred to other models. In conclusion we ought to comment on the formulation 
of a forward equation. We assume that the state of the system is known at t = 0, and our aim 
through this model is to predict the possible future states of the system and provide information on 
their statistical properties. This problem will be a useful model for a number of physical situations 
where the Burgers equation arises. 

We now introduce the backward problem. Here we assume that we know the final state of the 
system, which is the square integrable J^-measurable random field p : M. x 0, — >■ M, and we wish 
to find which initial condition must be chosen in order to drive the system at time T to this state. 
The proper formulation of this problem is in terms of the BSPDE: 



dU 



1 



a 2 (t, x)U xx + a(t, x)UU x + g(t, x)U^ u + b(t, x)U x 



dt 



(5) 



+ e{t, x)U + s(f, x)^> x + m(t, x)^ u + f{t, x) 

+ ^ u dW(t), 0<t<T, i£R, 
U(T,x) =p(x), iel. 

It should be stressed here that now the unknowns are the pair of random fields (U, and 
not just U. Therefore in contrast to what holds for the forward problem ([3]), ty(t,U) can no 
longer be thought of as given data to the problem, i.e. as pre-described by the model, but now it 
has to be specified as part of the solution of the problem. What is even more interesting is that 
ty u is uniquely determined by U, as a properly defined functional of U and not necessarily as a 
function of U. In some sense ty u has to be interpreted as the (unique) auxiliary process needed 
to drive the system to the desired random final state. Therefore, it may have the interpretation 
of a control procedure which acts on the system so as to drive it to a desired state. Furthermore, 
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the compensation of this auxiliary process, allows us to obtain a solution to the backward problem 
which is adapted to the filtration generated by the Wiener process. 

The important motivation for this problem arises from its connection with control problems. 
By the discussion above, one observes immediately that the backward problem is intimately related 
to the problem of controllability of the stochastic Burgers equation. This problem can be stated 
as follows: given an initial state can we drive the system by the proper control procedure to a 
desired final state which is random but with prescribed statistical properties? As we have stated the 
problem, it allows us to characterize this initial state that will drive us to the desired final state, and 
at the same time through ty u characterize the control procedure needed in the accomplishment 
of this task. The connection is much deeper; in particular, through the generalization of the 
Pontryagin maximum principle in the context of SPDEs one may show that the dual (adjoint) 
system associated with a large variety of optimal control problems will have a form closely related 
to ©. 

2.2 Notation, functional setting, and some important preliminary results 

Throughout this paper, denote by C fc (M n ) for any integer k > the set of real-valued functions on 
M n that are continuously differentiable up to order k. In addition, for any 1 < p < oo, any Banach 
space X with norm || • ||x, and any sub-u-algebra Q C J 7 , denote by 

• Lg($l,X) the set of all X-valued, ^-measurable random variables X such that < oo ; 

• Lp(0,T;X) the set of all F-progressively measurable, X-valued processes X : [0,T] x Q — >• X 
such that f q \\X(t)\\^dt < oo, a.s.; 

• Lp(0, T; L p (0; X)) the set of all F-progressively measurable, X-valued processes X : [0, T] x 
n ->■ X such that J Q T E\\X(t)\\^dt < oo ; 

• Cf([0,T];X) the set of all continuous, F-adapted processes X(-,u) : [0,T] — > X for P-a.e. 

Moreover, define similarly the set Cf([0, T]; L p (f2; X)) , and let M + stand for the positive real num- 
bers. Finally, for every random field F : [0,r]xR n xfi^M denote by A F , ^ F : [0, T) x R n x Q -> R 
the random fields for which F obtains the semimartingale decomposition 

F(t,x) =F(0,x) + / A F (s,x)ds+ [ ^ F (s,x)dW(s) V (t, x) € [0, T] x R n , 
Jo Jo 

almost surely. Assume a similar notation for stochastic processes as well.. 

In order to carry out computations regarding the change-of-variable formula for random fields, 
we recall the following useful implication of the generalized Ito-Kunita-Wentzell (IKW) formula 
(e.g. [31], Section 3.3, pp 92-93). 
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Proposition 2.1. Suppose that the random field F : [0, T] xR n x SI -> R is o/c/ass C°> 2 ([0,T] xR n ) 
and A F ,^ F : [0,T] x R n x O ->■ R are ^-adapted random fields of class C°' l ([0,T] x 1") and 
C°' 2 ([0,T] x R n ), respectively. Furthermore, let X = (XW, . . . , XO 1 ))* be a vector of contin- 
uous semimartingales, where A x 1 (•) is an almost surely integrable process and is an 
¥ -progressively measurable, almost surely square integrable process. Then F(-,X(-)) is also a con- 
tinuous semimartingale, with decomposition 

F(t,X(t)) =F(0,X(0)) + f A F (s,X(s))ds+ f* * F (s,X(s))dW(s) 

Jo Jo 

+ £/, ^rAs,X( s ))A^\s)ds 

+ £ f ^-F(s,X( S ))^\s)dW(s) 
i=i Jo 1 

+ E/ £-* F (s,X(s))* x{ '\ S )ds (6) 

i=i j=i j 

for every < t < T . 



3 Connection of Backward Burgers Equation with FBSDEs and 
a Probabilistic Approach to the Cole-Hopf Transformation 

For a given function p : R —¥ R and a ^ 0, consider the backward deterministic Burger's equation 
of the form 

U t + \o- 2 U xx -a 2 UU x = 0, 0<t<T, xGR, 

I \J) 

U(T,x)=p(x), xGR. 



It is well known that using a transformation to a new variable V : [0, T] x R — > R + such that 

dx 

1 

2 7 (8) 

^(r, x) = e" ^ (x)cte = q(x), x € R. 



[/ = — ^ lnT^, equation ([7]) is linearized assuming the form of the heat equation 

Vt + z o- 2 V xx = 0, < t < T, x G R, 



This is the celebrated Cole-Hopf transformation, through which one may construct solutions of 
Burgers equation using appropriate solutions of the heat equation. 
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Let us revisit this linearization from a probabilistic viewpoint. For any x € R, consider the 
system of FBSDEs 

dX(t) = adW(t), 0<t<T, 

dY(t) = aY(t)Z(t)dt + Z(t)dW(t), < t < T, (9) 
X(0) = x, Y(T)=p(X(T)), 

where X(-) is the forward process, Y(-) is the backward process, and Z(-) is the auxiliary process 
needed for the well-posedness of the problem, each defined on [0, T] x Q. If one looks for a 
Markovian solution to this problem of the form Y(t) = U(t,X(t)) then a simple application of 
Ito's formula yields both that the deterministic function U should satisfy equation ([7|) and the 
relationship Z(t) = aU x (t,X(t)). Consider further the FBSDE 

dx(t) = adW(t), 0<t<T, 

dy(t) = z{t) dW(t), < t < T, (10) 
x(0)=x, y(T)=q(x(T)) 

for the F-adapted processes x,y,z : [0,T] x — >• R, and follow the same steps as above to verify 
that for a Markovian solution of the form y(t) = V(t,x(t)) the deterministic function V should 
satisfy equation (jHJ) and z[t) = aV x (t,x(t)). 

The following issue is addressed. Taking into account that the PDEs (|7|) and ([Sj) are related 
via the Cole-Hopf transformation, we should expect that the probabilistic systems Q and (|10p . 
giving rise to these equations respectively, must be related as well. Conversely, if we find a point 
transformation between the variables (X, Y, Z) and (x, y, z) then we could use it to reproduce the 
Cole-Hopf transformation that allows us to get from the Burgers equation to the heat equation. 
In particular, we consider the point transformation between the solutions of ([9]) and (|10p : 

X = x, 

Y = y(x,y,z), (11) 
Z = Z(x,y,z), 

where y,Z : R 3 — > R are deterministic functions to be determined. Then we have the following 
result. 

Proposition 3.1. Suppose that the process triplet (x,y,z) solves U0\) . and the function y belongs 
to C 2 ' 2 ' 2 (R 3 ) and satisfies the PDE 

2 0-2 + ^ z 2 y yy + - h 2 y zz + a z y xy + hzy yz + ah y xz (12) 

-o- 2 yy x -o- z yy y -o-hyy z = o 

for every h € R. Furthermore let 

Z(x(t),y(t),z(t)) =o-y x (x{t),y(t),z(t)) +z(t)y y (x(t),y(t),z(t)) (13) 
+ h(t)y z (x(t),y(t),z(t)), 



10 



in terms of the process h : [0, T] X Q — )■ R suc/i that z(-) has the local martingale representation 
dz{t) = h(t) dW(t); cf. Remark \3.SX Then the process triplet (X,Y,Z) defined through the point 
transformation of Ul\) solves ([Pp. 

Proof. Take the Ito differential of Y(t) = y(x(t),y(t),z(t)) to get 

dY(t) = I(x(t), y(t), z(t); h(t)) dt + K(x(t),y(t), z(t); h(t)) dW(t), (14) 



where 



111 

I(x, y, z;h) = - a 2 y xx (x, y,z) + - z 2 y yy (x, y,z) + - h 2 y zz (x, y, z) 

+ o-zy xy (x, y, z) + hzy zy (x, y, z) + ahy xz (x, y, z), 
K(x,y,Z]h) = ay x (x,y,z) + zy y (x,y,z) + hy z (x,y,z). 



Thanks to (|12p and (|13p . these functions satisfy the relationships 

I{x{t),y(t),z{t); h(t)) = ay(x{t),y(t), z(t)) Z(x(t), y(t), z(t)) , 
K(x(t),y(t),z(ty, h(t)) = Z(x(t),y(t), z(t)) , 

and substituting them to (fl~4"|) we obtain the stated result. □ 

Remark 3.2. Note that applying Ito's differential to z(t) = aV x (t, x(t)) , in conjunction with 
allows us to take dz(t) = h(t) dW(t) in the previous proposition. 

Remark 3.3. One can easily verify that a solution of equation MSfy for any h G R is the rational 
function y(x,y,z) = —z/ay. This interpreted in terms of the solutions U and V of ([?]) and 
respectively, corresponds to the Cole-Hopf transformation, since 

u( t , * W ) = y W = y (,(<), m, m = = • 

Hence, making use of probabilistic tools we managed to derive an alternative representation of this 
transformation in terms of a solution of the PDE 



4 Linearization of the Forward Random Burgers Equation 

In this section we consider the forward version of the stochastic Burgers equation using the 
notation ^(t, U) = ^ u and stress the fact that here \P is a pre-described function of U. It is the 
object of the subsequent result to find the general form of the random coefficients of the above 
equation as well as of the function ty u such that this forward SPDE (FSPDE) can be transformed 
to a linear heat equation through the use of the Cole-Hopf transformation. Therefore, if we obtain 
somehow positive solutions for the linear one we may obtain solutions for the original Burgers 
equation with random coefficients by means of the Cole-Hopf transformation. 
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Theorem 4.1. The FSPDE 

1 



dU = 



a 2 (t,x)U xx + a(t,x)UU x + g(t,x)U^ u + b(t,x)U x + e(t,x)U 



+ s(t, x)^ u x + m(t, x)V u + f(t, x) 



dt 



(15) 



+ y u dW{t), 0<t<T, iGR, 
[7(0, x) =p(x), x€R 

is linearizable through the Cole-Hopf transformation to a stochastic heat equation if a, a, g, and s 
are stochastic processes of the form 



a (-) - °" 2 (-)> 9(-)-o-{-), and s(-) = -cr(-), 

t/ie function ty u is given by 

y u (t, x ) = £(t)U x {t, x) V (t, x) € (0, T]xt 

/or a stochastic process I: [0, T] x £7 — > R defined by 

<(-) = -2a(-), 



and finally 



e(t, x) = b x (t, x) + £(t)m x (t, x) V (i, x) G [0, T] x R ; 



in particular, if U G C F ([0, T]; L 2 (S1; C 2 (R))) zs a random field that solves FSPDE 



(16) 



(17) 



(18) 



(19) 



| a 2 ^)^ - a\t)UU x + ( 6(t,x) - 2a(t)m(t, x) ) 17 X 



+ e(t,x)E/ + /(t,x) 
[7(0, x) =p(x), xGR, 



+ £(t)£/ x oW(t), < t < T, x G R, (20) 



£/ien £/te random field V which is given by the Cole-Hopf transformation is of class 
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C F ([0,T];L 2 (ft;C 3 (R;IR + ))) and solves FSPDE 



<n <i ? cj 2 (t)T4.x + fc(<> a?)^„ - c(i, x)y }> dt 



+ £(t)^dW(<)i < t < T, x G 



(21) 



F(0,X) = e -/P(^ 4g( S ), X 



where 



c(t,x) = / f(t,x)dx + c(t) and k(t,x) = b(t,x) + £(t)m(t,x) 



(22) 



/or every (t,x) G [0,T] x R and some process c(-). Furthermore, the latter admits a solution of the 
form 

V{t,x) = G(t,x + H{t)), 0<t<T, xeR, (23) 



where 



H(t) = / ^(s)dW(s), 
'o 



< t < T, 

in terms of a positive random field G of class Cf([0, T]; L 2 (fi; C 3 (M; R + ))) that solves FSPDE 

dG =\\ o- 2 (t)G xx (t,x) + k(t,x - H(t))G x (t,x) 

- c(t, x - H(t)) G(t, x) | dt on (0, T] x R, 
G(0,z)=g(z) on R. 



(24) 



(25) 



Proof. We first establish the definitions of (|16p -(|19 p . according to which the FSPDE (|20p trans- 
forms to the FSPDE (|2ip . To this end, consider the Cole-Hopf transformation U = — J^lnV, 
in terms of a strictly positive random field V G Cf([0, T]; h 2 (Q; C 3 (R; R + ))) with a semimartin- 
gale decomposition such that ^ v G Lf(0, T; C 3 (R)). For every < t < T, Ito's formula and the 
uniqueness of the semimartingale decomposition imply that 



dU 



d_ 

dx 



dV (* 



V\2 



V 2V 2 



-dt 



and * 



V 



dx \~V~ 



(26) 
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and substituting back into (jT5j) we have 



d_ 

dx 



dV 



V\2 



2V 2 



-dt 



Id 3 d d 2 

-a\t,x)— 3 lnV + a(t,x)-lnV— 2 lnV 



. , d , T a f^ v 
^ x) dx l * V dxW 



d 2 



d 

e(t, x)—- In V 
ox 



m(t, x] 



d_ 

dx 



+ f(t,x 




(27) 



subject to the initial condition of (|2ip . Inspired by the methodology implied by the Cole-Hopf 
transformation in the linearization of ((TJ), we accordingly integrate the last equation with respect 
to x and assume that a, a, g, and s are independent of the spatial variable, i.e. are stochastic 
processes, to arrive at 



dV 



<9 2 /" V x 

b(t,x)—-^lnVdx+ / eft, x)— dx 
ox J V 

+ / m(t ' x) ^(V) d \ 



(28) 



2y 



This is an FSPDE for the random field V where ty v is a function of V, the exact form of which 
depends on the choice of the function ^> u through the second equation of (|26p . This is in general 
a nonlinear FSPDE, whose form can be as complicated or even worse than the original Burgers 
equation. However, to end up with a linear FSPDE (by analogy to what happened in the deter- 
ministic case where it reduced to the heat equation fl2J)), it suffices to consider the special ansatz 
of (|17p and Q18|) for the function VP , which also gives that 



x) = £(t)V x (t, x) V (t, x) e (0, T] x 



(29) 



We show now that this choice of the function ^ u (leading to the equivalent choice (|29jl of the 
function \P ) is enough to turn FSPDE (I28p for the evolution of V into a linear one of the heat 
type. Indeed, employ equation ([29|) to reformulate the integral expressions of (f28|) as 



y dx 



—— — — — — dx 



v 



2 
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and, in conjunction with integration-by-parts, as 



J Kt^)-^\nVdx + J e(t,x)— dx + J m{t,x)— I — I 



dx 



= b(t,x) + £(t)m(t,x) 

Hence, FSPDE (ggj becomes 

' 1 



Vx 

V 



+ 



e(t, x) - b x (t, x) - £(t)m x (t, x) 



V 



dx. 



dV= )2 O 7 ^ + 2s ® e ®) Vxx + ^ + ^) m (*> x 



Vr 



l r 



2 



(T 



; (t) + a(t) + ( 5 (t) + 2 S (t))^)-£ 2 (t) 



v 2 
v 



-dx 



e(t, x) - b x (t, x) - £(t)m x (t, x) 

-V J f(t,x)dxjdt + £(t)V x dW(t). 

In order for this equation to be linear with respect to the unknown field V, the coefficients in front 
of the terms V 2 jV and V x /V should vanish. Therefore, we consider that the stochastic process 
coefficients of FSPDE (|15p are expressed in terms of a(-) and, in particular, are defined as in (|16H . 
Regarding the coefficient of the term V x /V, this is the case by assuming either £{■) = a(-) or 
£{■) = — 2<r(-). We reject the first choice, which leads to definition (118p . since it reduces the last 
FSPDE into an ill-posed parabolic problem; specifically, for this choice the parabolic term comes 
with a negative diffusion coefficient. We may further assume the definition of (|19p to eliminate the 
coefficient of the term V x /V and finally reach the forward linear stochastic heat type equation of 
(I2T1) . given the definitions of (1221) . 

To prove the second assertion of the theorem, we seek for positive solutions of (|2ip that have 
the form of ([23]). Here, we assume that G is a random field of class C F ([0, T); L 2 (ft; C 3 (M; R+))) 
with a semimartingale decomposition in terms of the pair (A G ,^ G ) £ Cjf([0, T];L 2 (il; C 1 (1R))) x 
L|(0, T; C 2 (M)) . Employing the generalized IKW formula to (|23p . in conjunction with (|24p . we 
obtain 

dV = ^A G (t,x + H(t)) + l - G xx [t, x + H(t)) (\P H (i)) 2 

+ V H (t)V G (t,x + H{t))\dt 
+ {* H (t)G a (t, x + H(t)) + * G (t, x + H(t))}dW(t) 
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for every (t, x) € (0, T] x R. Making also use of (f23j) . we rewrite (f2Tj) as 

= j|(7 2 (t)G !M (t,:E + .ff(t)) +k(t,x)G x (t,x + H(t)) 

- c(t, x)G(t, x + H{t))\dt + l{t)G x (t, x + H(t))dW(t), 
for every (t, x) € (0, T] x 1. A comparison between the last two equations leads to the definitions 

A G (t,x) = - a 2 {t)G xx (t,x) + k(t,x - H(t))G x (t,x) - c(t,x - H(t))G(t,x), 
* G (v)=0, 

where (|18p has been used as well. Finally, substituting them in the semimartingale decomposition 
of G , we derive the FSPDE of fj25[> : the initial condition comes from (|2ip and (|23p for t = 0. 
Regarding the positivity of G see the following remark. □ 

Remark 4.2. FSPDE A25\) is a linear PDE with random coefficients, but does not contain any 
stochastic integral. This means that it can be treated pathwise rather than in an ltd integration 
sense, a fact that simplifies its analysis immensely, and allows us to use qualitative results, e.g. 
maximum principles, to study properties of its solutions like positivity, monotonicity etc. In fact, 
for the special case of constant coefficients with k(-, •) = c(-, •) = this PDE reduces to 



5 Linearization of the Backward Random Burgers Equation 

In Section [3] we showed how the backward deterministic Burgers equation can be linearized using 
its connection with a system of FBSDEs, through which we obtained an alternative probabilistic 
derivation of the celebrated Cole-Hopf transformation. In the present section we address the 
question of whether this analysis can be generalized for the backward stochastic Burgers equation 
with random coefficients of ([5|). It will turn out that in order to answer this problem it is no 
longer enough to use the Cole-Hopf transformation that works for the deterministic or the forward 
random problem but rather we have to extend it in such a fashion as to take into account the 
intricate nature of the backward problem, and the presence of the related unknown process pair 
(U,V U ). 

For the simplicity of our mathematical program which aims to the linearization of ([5]), we also 
assume that the coefficients are bounded and infinitely differentiable with respect to the spatial 
variable x, with all their partial derivatives bounded as well. Carrying out though our analysis, 
we shall see that both weaker regularity conditions in x may be imposed and explicit relationships 
among the coefficients will be in need. 

5.1 Connection with FBSDEs of random coefficients 

We consider the forward process X(-) that satisfies X(t) = x + £ a{s,X(s))dW{s) for (t,x) € 
[0,T] x R, and seek for the corresponding backward given by Y(t) = U(t,X(t)) in terms of a 
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solution (U,V U ) 6 C F ([0,T];L 2 (n;C 3 (M))) x L|(0, T; C 2 (M)) of BSPDE ©. Note here that 
we have considered an additional order of smoothness for the pair (U, than is required to 
constitute a solution of ([5]). Nevertheless, this extra assumption allows the application of the 
generalized IKW formula, which coupled with ([5]) yields 



dY(t) 



a(t, X(tj) U x (t, X(t)) + g(t, X(t)) * u [t, X(t))] U(t, X(t)) 

+ b[t, X(t)) U x (t, X(tj) + e(t, X(t))U(t, X(t)) 

+ (s(t,X(t)) +a{t,X(t))^(t,X(t)) +m{t,X(t)p u (t,X(t)) 

+ f(t, X(t)) Xdt + \a(t, X(t)) U x (t, X(t)) + V u (t, X(t))]dW(t), 



for < t < T. In order both to be consistent with the form of the deterministic Burgers equation 
of ([7]) and to avoid matters of technical fuss in our following analysis, we shall consider that 
g(-) = <t(-) is a stochastic process. Additionally, by setting 



-a{t) 



(30) 



a(t,x) = g(t)a(t), b(t,x) = a(t)m(t,x), and s(t, x) = 
for (t,x) 6 [0,T] x R, we obtain the first part of the following result. 

Proposition 5.1. (i) Suppose that the pair of random fields (U, is of class CV([0, T]; L 2 (fi; C 3 (R))) 
xL|(0, T; C 2 (M)) and solves the stochastic Burgers type BSPDE 



dU 



\ o 2 {t)U xx + a 2 (t)UU x + a(t)Um u + a(t)m(t, x)U x 



+ e(t, x)U - <r(t)*£ + m(t, x)^ u + f(t, x) 

+ $> u dW(t), 0<t<T, xgR, 
U(T,x)=p(x), i£l. 

Then the triplet (X, Y, Z) of stochastic processes, given by 



dt 



X{t)=x+ a(s)dW{s), Y(t) = U(t,X(t)), 
Jo 

and Z(t) = a{t)U x (t, X(t)) + <1> U (t, X{t)) 
for every t € [0, T], satisfies the FBSDEs 

dX(t) =a(t)dW{t), 0<t<T, 

dY(t) = \a{t)Y{t)Z{t) + e(t,X(t))Y(t) + m(t,X(t))Z(t) 

+ f(t,X(t))]dt + Z(t)dW(t), 0<t<T, 
X(0)=x, Y(T)=p(X(T)). 



(31) 



(32) 



(33) 
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(ii) Consider that the pair of random fields (V, \P ) belongs to the class 
C F ([0,T];L 2 (O;C 4 (]R))) x L|(0, T; C 3 (R)) and satisfies the stochastic heat type BSPDE 

dV = {~\ ^\t)V xx + a(t)d{t, x)V x - <r(i)tf£ + d(t, x)V v 

+ c(t,x)v}dt + ^ v dW(t), 0<t<T, ieR, (34) 
V{T,x) = q(x), iGR. 

Then the triplet (x, y, z) of stochastic processes, defined by 



(35) 



x(t)=x+ a(s)dW(s), y(t) = V(t,x(t)), 
Jo 

and z{t) = a{t)V x (t, x(t)) + ^ v (t, x(t)) 
for every t € [0, T], solves the FBSDEs 

dx(t) = a(t)dW(t), < t < T, 

dy(t)= c(t,x{t))y(t) + d(t,x(t))z(t) dt + z(t)dW{t), < t < T, (36) 
x(0)=x, y(T) = q(x(T)) . 

Proof. We need only to show (ii). Indeed, the BSPDE ([36]) follows readily through (j34"J) . (f35j) . and 
the generalized IKW formula. □ 



The first part of this proposition provides a solution to the FBSDEs (J33J) by means of a solution 
of the Burgers BSPDE (|3ip . In Section[3]though, we saw that the Cole-Hopf transformation reduces 
the Burgers equation (|7|) to the linear heat equation ([8]). Therefore, its next statement addresses 
the most general type of a backward linear stochastic heat equation (cf. (|34p ) that is associated 
with an FBSDE system (cf. (pBl ). 



5.2 Characterization of the generalized Cole-Hopf transformation in terms of 
an SPDE 

Our goal in this subsection is to investigate wether there exists a stochastic version of the Cole-Hopf 
that reduces (|3"Tj) to (f3"4"j) for appropriate random fields c,d : [0,T]xRxfl->R and q : R x $7 — > E. 
According to Remark 13. 3\ in the deterministic case the Cole-Hopf transformation is represented 
by the point transformation (|lip . Thus, by analogy we shall seek for a point transformation of the 
form 

X(t) = x(t), 

Y(t) = y(t,x(t),y(t),z(t)), (37) 
Z{t) = Z(t,x{t),y{t),z(t)), 
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that relates the solutions (X, Y, Z) and (x, y, z) of (f33|) and (f36|) respectively, in terms of suitable 
random fields y, Z : [0, T] x M 3 x £1 — > M. This transformation will provide explicitly the Cole-Hopf 
one in the generalized case of stochastic coefficients. 

Assumption 5.2. Let the stochastic coefficient a (■) obtain a semimartingale decomposition, where 
A a {-) is an integrable a.s., ¥ -progressively measurable stochastic process, and is a square- 

integrable a.s., ¥ -progressively measurable one. Furthermore, let the latter process obtain a similar 
semimartingale decomposition in terms of the stochastic processes A^" (•) and (•), respectively. 



Proposition 5.3. Consider Assumption \5.2[ Let also (x,y,z) be a solution of FBSDEs \3b}) . 
(y,ty y ) be a random field pair that belongs to C F ([0, T]; L 2 (ft; C 3 (M 3 ))) x L|(0, T; C 2 (M 3 )) and 
satisfies the BSPDE 



dy 



2 °*m 



XX 



zz 



o-(t) z y 



o-(t)hy 



xz 



zhy 



+ o-\t)yy x + a{t) z yy y + a(t)hyy z + o-(t)m(t, x)y 



X 



(38) 




v 





(t,x(t)) 





on 



[0, T) x M 3 



for any h & R, and 



z±m y + a{t)y x + z(t)y y + h(t)y z 



(39) 
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evaluated at (t, x(t),y(t), -z(i)) on [0, T). Then the process triplet (X, Y, Z) given by the point trans- 
formation £Ffy solves the FBSDEs 

Given 



Proof. Firstly we provide the semimartingale decomposition of the process z{-) in 
Assumption 15,21 considering that the pair (V, belongs to the class Cf([0, T]; L 2 (J7; C 4 
Lp(0, T; C 3 (IR)), and assuming that the random field has a semimartingale decomposition 
for (A* V ,^ V ) e C F ([0,T];L 2 (n;C 1 (M))) x L|(0,T; C 2 (R)), employ the product rule and the 
generalized IKW formula to process z(-) of (f35|) . in conjunction with (j34l) . to get that 

dz{t) ={A* V (t,x(t)) + a(t)(o-(t)d(t,x(t)) + ^(t))v xx (t,x(t)) 



+ 



A a (t) + a 2 (t)d x (t, x(t)) + a(t)c(t, x(t)) 

1 



V^(t,x(t)) 



+ <x(i) Ccc (t, x(t))V{t, x(t)) + - a 2 (t)^ x (t, x(t)) 
+ (a(t)d(t,x(t)) +**(t))*Y(t,x(tj) 



+ a(t)d x (t,x(t))^ v (t,x(t)) +a(t)^ x V (t,x(t)) \dt 



+ 



^ V (t,x(t)) + a 2 (t)V xx (t,x(t)) 

+ ^{t)V x {t, x(t)) + 2*(t)*v (t, x(t)) 



dW{t) 



for < t < T. Therefore, set 

h(t)±* 9V (t,x(t)) +a 2 (t)V xx (t,x(t)) + V a (t)V x (t,x{t)) + 2a(t)^(t,x{t)) 
for every t £ [0, T] and recall the third equation of (|35j) to derive the relationships 

z(t)-* v (t,x(t)) 



V x {t,x(t)) 
V xx (t,x(t)) 



a{t) 



H t)-^ v (t,x(t)) * g W(j* v M'))-*(') 



a 2 (t) 



a{t) 



(40) 
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and deduce eventually for any t £ [0, T) the semimartingale decomposition 



dz(t) 



A a (t)+a 2 (t)d x (t,x(t)) - $> a (t)d(t,x(t)) 



(**(t)) s 



a(t)c(t,x(t)) 



£(*) 
a(t) 



+ \o 2 m v xx {t,x(t)) - (a(t)d{t,x(t)) + tt ff (i)W(t, a? (t)) 



A a {t) - ^ a (t)d(t,x(t)) 
+ a(t)c(t,x(t)) - 



(41) 



a(t) 
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a{t) 



+ a* v (*) + (7(t)*r («, »w) - ( ^w) + ^32 ^ 



+ <r(t)ci B (t,x(t))y(t) + 



a(t) 



+ d(t,x(t)) )h(t) }dt + h(t)dW(t). 



We are ready now to argue that BSPDE (|38p constitutes a sufficient condition for the point 
transformation (|37|) to produce solutions of the FBSDEs (j33l) directly out of solutions of the 
FBSDEs (|36p . Indeed, this follows directly from the generalized IKW formula applied to the 
second equation of (|37l) . in combination with (j36j) . (JHJ), and (|38l) , □ 



We shall focus next on establishing an equivalent characterization for the solutions of BSPDE 
(I38p . in order to determine the point transformation of (|37j) . A byproduct of this analysis will 
be the Cole-Hopf transformation between the solutions U and V of the BSPDEs (|3ip and ([34 
respectively; cf. Remark 15.51 



(42) 



Proposition 5.4. Lei Assumption \5.2\ hold and set 

d(t,x) = m(t,x) on [0,T] x M. 
Then the BSPDE of A38\) admits solutions of the form 

y(t,x,y,z) 
on [0,T) x R 3 , where 



a(t)y-r(t,x) a 2 {t) 



W{t) a(t)r x (t,x)-^(t)y + W(t,x) 

T~ o / . \ "T" 



a(t)(a(t)y-r{t,x)) 



(43) 



e(t,x) = a(t)m x (t,x) and c(t,x) = — J f(t,x)dx + c(t) on [0,T] 



(44) 
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for some stochastic process c(-), the random field r : [0, T] x R x -> R satisfies the linear BSPDE 
dr=\ - i a 2 (t)r xx + (v 7 (t) + o-(t)m(t,x)y x 



+ 



+ 



a{t)c(t, x) - y a {t)m(t, x) - 2 + A a {t) 



a{t) 



a(t) 



a(tW x (45) 



^ a (t) 



m(t,x))y r }dt + y r dW(t) on [0,T)xK, 



and the constraint 



2(^(t)Y , ~ 



a(t) 



+ (t) 



2(^(t)) 2 m(t,x) 2^(t)^(t) ^ " 

+ W) + () 



+ 



a(t) a 2 (t) a{t) 



a(t) 



a(t) 



-a(t)*f + + m(t, x)) + I <r 3 (i)* 

-o^)(a(i)m(t,x) + * CT (i)W 



v 

XX 



A a (t) - ^ a (t)m{t, x) + a(t)c{t, x) 



+a(t)A* V + <7 2 (t)** V - ( a(t)m(t, x) + ^ V - A w = 



(46) 



hobb on [0, T) x R. 

Proof. Since PS|) holds for all h G R, we see that from the term quadratic in /i we get 3^2 = 0, 
which on the real line shows that 



y(t,x,y,z) =V(t,x,y)z + Q(t,x,y) on [0,T)x 



(47) 



for appropriate random fields P,Q6 Cf([0, T]; L 2 (Q; C 3 (]R 2 ))) , that each obtains a semimartingale 
decomposition. Additionally, from the terms linear in h and ()42j) we find that 



y z -^ = 0; 
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substituting (|47p . the left-hand-side becomes a linear polynomial with respect to z which implies 
that 



P(t,x,y) 



1 



a{t)y - r(t,x) 



and 



(48) 



sd{t, x, y) — n , - + 



a(t){a(t)y-r(t,x)) 



(49) 



on [0, T)xl 2 , in terms of a random field r : [0, T] xlx Q — > R with a semimartingale decomposition 
to be determined. 

As a consequence of this analysis, the BSPDE (j38|) takes the reduced form 



dy 



A a (t) + a 2 (t)m x (t,x(t)) 

-^ a (t)m(t,x(t)) + a(t)c(t,x(t)) 



(W(t)) 2 \z(t) 



a(t) Ja(t) 



+ 2 " ( <r(i)m(i,z(i)) + tt'(t) ) ^ (t, x(t)) 



A a (t) -^ a (t)m(t,x(t)) 

+ ^)c(M(*))-^^*M) +A .^) 



a(t) 



er(t) 



(50) 



+ a(t)^* V (t,x(t)) - (m(t,x(t)) + 



a(t) 

+ a(t)c x (t,x{t))y(t) 



y z 



+ a (t)m(t, - c(t, x)yy y + e(i, x)y + <r(t);ytf y - cr(t)*^ 
- z^ y + m(t,x)^ y + f(t,x) \dt + V y dW(t) on [0,T)xR 3 



and from (|4"7|) . (jl5]l . and (f4Uj) obtains solutions in the form of ([4~3|h Through substitution of (|4"7|) 
and comparing the finite- variation terms of the two sides, we arrive again at an equation between 
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two linear polynomials with respect to z which yields the relationships: 



A v 



- a(t)m(t, x)V x - c(t, x)yV y + e(t, x)V + a(t)V^ Q + a(t) Q® v 



<r 2 (t)V xx - a(t)Q xy + a 2 (t)QV x + a 2 (t)VQ x + a(t)QQ y 
V 

A a (t) + a 2 (t)m x (t, x) + a(t)c(t, x) 



^ a {t)m(t,x) 



(*'(t)) s 



a{t) 



V 

a(t) 



a(t)V% -V$ + m(t,x)V v , 



(51) 



A Q 



1 



v 2 (t) Qxx +o{t)[a{t)Q + m(t, x))Q x - c(t, x)y Q y + e(t, x) Q 



- a 2 (t)^ x (t,x(t)) - a(t)m(t,x(t)) +^(t)W(t,x(t)) 



A a (t) - $> a (t)m(t,x(t)) 
+ a(t)c(t,x{t)) - 



V T (t)) 2 \ ^ v (t,x(t)) 



a(t) 



a{t) 



+ A* V (t) + <r(t)*r (t, x(t)) - ( m(t, x(t)) + ) ^ V 



a(t) 



+ a(t)c x (t,x(t))y 



V - o{t)V% + (a{t)Q + m(t, x)\ * Q + f(t, x) 



both on [0, T) x M 2 . Employ ltd differentials to (l48j) and (f49|) to compute explicitly the semi- 
martingale decompositions of V and Q, and make the substitutions in (|5ip to end up at two 
zero-polynomials in y, of first and second order respectively, which imply eventually the definitions 
of flMJ), the BSPDE (jig]) , and the constraint (jS}. □ 



BSPDE (|45p admits the apparent solution of r(- , •) = 0; then, (I43D simplifies to 



y(t,x,y,z) 



<?{t)y 



on [0, T) x 



(52) 



as well as, solving (|46p for A* V and substituting back to the semimartingale decomposition of 
we get the BSPDE constraint 



d^ 1 



1 



cr 2 (^L+ a(t)m(t,x) + * ff (t) 



.r 



+ 



m(t, x) + c(i, x) 



a(t)* 



(53) 



dW(t). 
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on [0,T) x R. Moreover, combining (jglZj) . (|55) . ([57]> . and ([55]) we have that 



tf(t, x(t)) = r(t) = y(t, x(t), y(t), *(t)) = 



(54) 



V x (t,x(t)) * v {t,x(t)) 



V{t,x(t)) o{t)V(t,x{t)Y 



which is always the case if 



d ^ v 

U = -—lnV-—— (55) 
dx a(t)V y ' 

holds for every (t, x) € [0, T) x R. In fact for t = T, putting together the terminal conditions of 
(l3Tj) and (|34"1) . we get that the random field q must satisfy the condition 

P (x) = -— lng(x)- \ ' , xGR. (56) 

Remark 5.5. T/ie relationship of i55\) is the generalized Cole-Hopf transformation that provides 
the pattern to construct solutions (U,^ u ) of the stochastic Burgers BSPDE \31\) from solutions 
(V, \P ) of the stochastic linear heat equation subject to the additional linear constraint of 

\53\) and the terminal condition of t56\) ; a verification result is provided below (cf. Theorem 15. 6\) . 
Apparently in a deterministic setting \55\) reduces to the standard Cole-Hopf transformation of 
Section and A53\) becomes a tautology. 

One can immediately see that considering either deterministic or random coefficients the point 
transformations y of 177]) and (37\ ) are both of the same form; cf. Remark \3.3\ and [hlfy . However, 
the application of the probabilistic characterization of the Cole-Hopf transformation, deployed in 
the previous section, led to Proposition 15.31 and eventually concluded to the additional constraint 
for the random field ^> v of It 53]) . which vanishes in the deterministic case. 

Theorem 5.6. Consider Assumption \5.Sl definitions and |^[ ), and a random field pair 

(V,V V ) e Cf([0,T];L 2 (O;C 4 (R))) x L 2 (0, T; C 3 (R)) that satisfies the heat-type BSPDE Qlfl with 
V(-,-) > 0, subject to the constraints (S3D for the pair € C F ([0,T];L 2 (Q; C 3 (R))) 

xL|(0, T; C 2 (R)) and jl56\) for the random field q. Then the random field pair (U, 1 >i> u ) of class 
C F ([0,T];L 2 (S1;C 3 (R))) x L§(0, T; C 2 (R)) defined through {55P is a solution of the Burgers-type 

bspde jm 



Remark 5.7. It is important at this point to make a comment on the constraints as expressed in 
equations 153]) and 156\) . For a given final condition q(x), the backward heat equation \3J$ admits 
a unique solution pair (V, Then, the uniqueness of the semimartingale decomposition for 

the random field ^1> V provides us with a uniquely determined . Therefore, a fruitful way of 
handling these constraints is to substitute in them the pair (ty v , ), as defined by the solution of 
\34\ ), arid then find the general class of coefficients of the original problem so that the constraints 
are satisfied. This approach provides us with large families of exactly solvable backward Burgers 
equations with random coefficients; see Subsection \5.3l 
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Remark 5.8. It is worth mentioning that, regarding fijty) , j44)) > an d solution pairs (V, ) and 
(r, \& r ) of the BSPDEs (3~4\ ) and fijjfy, respectively, the random field y of Jff7| ) defined through 
provides the most general form of solution to [3$\) , subject to the constraint In view of (54] ) 

and h55]) we obtain here the relationship 

u = V x _ * v + <r(t)r x - *°(t)V + W 



V-r a{t)V-r a 2 {t) a(t)(a{t)V - r) 

for every (t, x) € [0, T) x R, which establishes the most general formulation of the Cole-Hopf 
transformation between the BSPDEs of \31\) and \3J$ - 

Proof. OF theorem ESI For < t < T and x G R, Ito's rule, ([55]), ([SD, and ([53]) yield 
dU = I a{t)^- - (a(t)m(t, x) + W(t)) + W(t)m(t, *)^jv 

+ ^-m(t,x)——-a(t)—-^+m(t,x" 



V y ' J a{t)V w 2V 2 y,J V 2 

(ibV)2 yV^V (^V\3 (i$ ,V)2 yVy* v 

+m ^k^-^-^-^ {t) hk- + ^ (57) 

2tj_\^xxx 2u\^ xx ^ x U\ 14. \^ x u\ 14. \^xx 

+ a(t)m(t, x)——^ m^i, x)—— — m(t, x)—— + m(t, x) 



+ o{t)-^ - a{t )^L - Cx (t, x) + _ \ dt 



Setting 



* - — r~\ — h , 9 +^ ,fT (t)^-- - + 



in conjunction with (155 1) . it is straightforward to obtain that 



^ v ^ V V ^ v ^ v (ty v ) 2 V ^ v 



a(t)V a( r )V 2 a(i)V 2 <7(t)V 3 v ' a 2 {t)V 



V 2 (t)F 2 v v 2 v 2 v 3 



a(t)V <j{t)V 2 V V 2 
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(V V ) 2 V X 


a 2 (t)V 2 


a 2 (t)V 3 




{V x f 


+ y2 


V 3 ' 




(^)3 


a 2 {t)V 2 


a 2 (t)V 3 



(tyV\2 ^V^V ^ V V 



_ A^yv* , ^ * v v* , _ 04-) 2 * v 

cr(t)y 3 l V(t)F 2 y 2 y 3 ' 

TT - _ l _ , _ o (^) 3 

a(t)y + a(i)F 2 V V 2 V 3 

a(t)1/ 2 C7(t)F 3 ' 

Therefore, recalling (|%2"]l and direct substitution of the above into ([3T|) leads to ([57|) and 

completes the proof. □ 

Theorem 15.61 necessitates the investigation of existence, uniqueness, and positivity of solutions 
for the linear BSPDE fj34f) subject to the constraints (153j) . which is another linear BSPDE for 
their martingale part, and (|56p . Ma & Yong ([31], [26]) studied linear BSPDEs of parabolic type 
and provided regularity conditions on their coefficients that ensured existence, uniqueness, and 
comparison results between solutions. However, their results are invalid under the presence of a 
BSPDE constraint. 

Corollary 5.9. Let Assumption \5.^ hold and V belongs to the class C^([0,T];~h 2 (Q; C 4 (M))) be a 
positive random field that satisfies the BSPDE 

dV = ^-^a 2 {t)V xx + a(t)m(t,x)V x + c(t,x)V^dt, < t < T, xeR, 
V(T,x) = q(x), iEK, 
where the random field c is given by an d Q * s defined as in \21\) . Then the random field 

U = -j^foV € C ¥ ([0,T];h 2 (n;C 3 
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defined through the Cole-Hopf transformation, solves the Burgers-type BSPDE 

1 



dU 



o- z (t)U xx + a z {t)UU x + a(t)m(t,x)U x + a(t)m x (t,x) 



+ a(t)m x (t,x)U + f(t,x) 



dt, < t < T, x € 



U(T, x) = p(x), x G 



Proof. The statement of the theorem follows as a direct application of Theorem 15.61 to the random 
field pair (V, 0); in particular, the constraint (|53|) holds trivially and (|56j) yields the definition of q 
as in (ED). □ 



5.3 Examples of completely solvable backward random Burgers equations 

In what follows, we present examples of random field pairs (V, ty v ) that satisfy the hypotheses of 
Theorem 15.61 for appropriate coefficients. 



Example 5.10. We postulate positive solutions for (|34|) of the form 

V(t,x)±ef 1 M w ® (58) 
on [0,T] x R, in terms of a given function f\ : R — > R of class C(R 4 ). ^From Ito's formula we have 

dV = h{x)VdW(t) + ^fl{x)Vdt, (59) 

which shows that 

^ v = h(x)V. (60) 
Then, the terminal condition of ([34p indicates that the random field q is determined by 

g(o;) = e /l(l)l1/(T) , x € R, (61) 



thus the random field p of (|31j) must be chosen according to (|56|) . Using (|58|) and (|60|) . we derive 
that 

K = /((x)viy(t), 
= f{(x)vw{t) + (/;(x))Vif 2 w, 

=/i(z)^ + /l(z)/i(z)VT^*), 
= f{{x)V + 2(/i(x)) Vw(t) + h{x)f'{(x)VW(t) 
+ f 1 (x){f[(x)) 2 VW 2 (t), 
and substituting them into (|34p . in comparison with (j59j) . we get 

\fl{x) = ~ \<r\t)[fl{x) + (f[(x)fW(t)] W(t) + a(t)m(t,x)f[(x)W(t) (62) 
+ m(t,x)/ 1 (x)-^)/((x)[l + / 1 (x)W(t)] + c(t,x). 
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Furthermore, Ito's formula and (J60J) give 



which implies that 



d* v = fl(x)VdW(t) + lff(x)Vdt, (63) 



= fl{x)V 

*r = 2f 1 (x)f[(x)V + /i(aO/i(x)VW(t); 



then, substitute these equations into ()53j) and compare it with (1631) to get 

\ff{x) = ~\ [f"W + 2(f[(x)) 2 W(t) 

+ /^/('(^(i) + h{x)(f[(x)) 2 W 2 {t) 

+ (<r(i)m(i,x) + /((a?) [1 + fi(x)W(t)] (64) 

/^(t) f± . . , (^ CT (t)) 2 \ „ , , 

- a(t)f 1 (x)f[(x)[2 + fi(x)W(t)] + ( K ^- + m(t,x) ) jf!(x). 

Obviously the random field triplet (V, * v , *® v ) belongs to the class C F ([0, T\; L 2 (ft; C 4 (IR))) x 
C F ([0, T]; L 2 (0; C 3 (M))) xL 2 (0, T; C 2 (M)) , according to TheoremESl Furthermore, equations UgSJ 
and (|64p form a 2 x 2 linear system which can be easily solved for the random fields m and c. A 
selection of these coefficients together with (fUTj) constitute the pair (V, Vl/ ) of (f58|) and (f60|) as a 
positive solution of (j34)) subject to the constraint (|53l) . 



Example 5.11. Alternatively, for a given positive function fa : R — > M + of class C(IR 4 ), we 
consider that (1341) admits a solution of the form 



V(t, x) 4 / 2 (x)e w/ W € C F ([0, T]; L 2 (ft; C 4 (M))) . (65) 

Now, Ito's rule yields 

dV = VdW{t) + - Vdt and ^ y = V, (66) 



the terminal condition of (|34p imposes that the random field g must be given by 

q(x) = f 2 (x)e w(T) , xeR, (67) 

and the constraint (|56p designates the random field p of (|31|) that should be considered. Therefore, 
differentiating ([65]) we obtain 

V x = &(x)e w ® and ^ SB = /^(s)e w W, 
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and through substitution back to (|34p and (|53p . a comparison with (|66p reveals that 

\h{*) = -\<J 2 {t)f' 2 \x) + a(t)m{t,x)f' 2 {x)+m{t,x)f 2 (x) 
-a(t)f' 2 (x) + c{t,x)f 2 {x), 

1 . , s 1 



2 



/2(x) = - ^\t)f 2 \x) + (a(t)m(t,x) + 

- <T{t)&{x) + + m(t, x)\ h{x) 

should hold respectively. Similarly to (i), the solution of this 2x2 linear system determines the 
proper random fields m and c for which the pair (V,tf v ) of <|E3) and (JBBJ) is a positive solution of 
subject to the constraint (f53l) . 



6 Stochastic Feynman-Kac Formulae 



In this section we establish stochastic Feynman-Kac results for the solutions of the Burgers FSPDE 
(|20p and BSPDE (|3ip ; cf. Propositions I6.ll and 16.21 below, respectively. Both results rely on the 
Cole-Hopf transformation, in its simple or generalized version (cf. (|55p) respectively, which reduces 
the initial SPDEs to linear ones. 



Proposition 6.1. Let the assumptions of Theorem \4-l\ hold, let the processes a, c and the random 
fields b, m be constants, let /(•,•) = 0, as well as p £ C 2 (R). //, in addition, there exist constants 
K > and < a < 1/2T such that 

max \G(t,x)\ + max \GM,x)\ < Ke ax2 Vi£l, 

0<t<T 0<t<T 

then the solution random field U of the corresponding Burgers FSPDE jl20\) admits for every 
(t,x) G [0, T] x R and a.s. the stochastic representation 



where 



G(t,x) 



U(t,x) 



q(x(t)) 



G x (t,x + H(t)) 
' G(t,x + H(t)) 



and 



G x (t, x) 



-Ct JjJX 



<lx{x(t)) 



here, the notation E x stands for the expectation that corresponds to the process 

dx(t) = kdt + adW(t), < t < T, 
x(0) = x. 
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Proof. According to hypothesis, differentiation of (|25p implies that the random field G x solves 
the same FSPDE with G but with initial condition given by G x (0,x) = q x (x) for x G R. Then 
the conclusion of the proposition is obtained by the Cole-Hopf transformation, (j23j) . and standard 
Feynman-Kac results (e.g. [28j Corollary 4.4.5) applied to the FSPDEs whose solutions are the 
random fields G and G T . □ 



Proposition 6.2. Consider the FBSDEs of \36\) on the interval [t, T] with x(t) = x for some 
< t < T, and assume that the hypotheses of Theorem \5.6\ hold with m{ • , •) =<£(•,•) = ^ a {-) = 0. 
Then the solution (U,^! u ) of the Burgers BSPDE \31\). corresponding to Ii55\) . obtains for all 
(t,x) S [0,T] x M and a.s. the representation 



U(t,x) 



z x (t) 
o-(t)y*(t) 



where 



z x (t) = E x 



+ 



-/ 1 T (^ 1 +4. I W)) rf V(T)p(,(r)) g (x(T)) 



^ + < 6 > x ^) de a( S )f(s,x(s))y*( S )d. ; 



y x (t) = E x \q(x{T))e-£< s ^ s))ds F{t) 



here, we denote by E x the expectation that corresponds to the forward process of the aforementioned 
FBSDE system. 

Proof. Thanks to hypothesis and the constraint (j53h . the BSDEs of (|36|) and (I4ip simplify to 

dy{s) = c(s, x(s))y(s)ds + z(s)dW(s), t < s < T, and 



dz(s) 



a(s) 



+ c(s,x(s)) I z(s) - a(s)f(s,x(s))y{s) 



ds + h(s)dW(s), 



for t < s < T, respectively. Denoting their solutions by y x (-) and z x (-), observe that these BSDEs 
are linear with drift terms independent of their auxiliary processes z(-) and h(-), respectively. 
Additionally, set t = T in (|54p and invoke the terminal conditions of (|3ip and (|36p to verify that 
the terminal condition of the last BSDE above is given by 



z(T) = -a(T)p(x(T))q(x(T)). 



(68) 



Therefore, the assertion of the proposition follows as an application of Corollary 6.2 by Ma & 
Yong (1997), in combination with (|54p . the above BSDEs, and their terminal conditions in (|36p 
and (1681). □ 
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7 Applications 



The following applications illustrate the use of Theorem 15.61 to construct solutions of backward 
stochastic Burgers equations with random coefficients, which are less commonly found in the liter- 
ature. These equations turn out to be associated with problems on controllability or mathematical 
finance. 

7.1 Controllability of a backward random Burgers equation 

Consider the following control system 



where (U\, U2) is a pair of control processes to be determined and po(x) is the initial state of the 
system. The question we wish to answer is the following. Setting a desired final condition p(x) for 
the control system (|69p . for which initial state po(x) does there exist a pair of control processes 
(Ui, U 2) such that this system can be driven in time T to the desired state p{x)l 

The linearization results for the backward Burgers equation may prove useful to answering 
this question. Assume that the two control procedures are connected by U% = -§^U\. Then the 
control system is equivalent to the BSPDE ([5]) with U\ = *f> u . Therefore, having obtained a 
solution (U, \P ) to problem ([5]) through the analysis of Section [5j we immediately obtain the 
control procedures needed to drive the system to the desired final state as U\ = ^ u and U 2 = ^ x , 
and we characterize the initial state we need to start from by po(x) = U(0,x). 

7.2 Pricing a contingent claim 

We consider a financial market on a finite time-horizon [0, T] consisting of a money market and a 
stock, whose prices Sq(-) and S(-) evolve according to the SDEs 



respectively. The interest rate r(-) > 0, the instantaneous rate of return //(•), and the volatility 
er(-) are taken to be bounded, F-progressively measurable random processes with a(-) ^ 0. 

In this market, we consider an economic agent with an initial endowment Y(0), who at any 
time t E [0, T] invests a proportion ir(t) of his wealth Y(t) in the stock, saves the remaining amount 
[1 — ir(t)]Y(t) in the money market, consumes with a predetermined rate C(t) > 0, and is obliged 
to pay taxes with rate L(t) > 0. Here the investor is not allowed to sell short stocks or to borrow 
money at the bond interest rate r(-); namely, the ir(-) above is restricted to take values only in 
[0, 1]. Furthermore, the tax regulation of the market mandates payments proportional to the size 
of the wealth Y(-) with respect to the interest rate r(-), but also applies in favor of investing to the 



dU 



- - a 2 (t, x)U xx + a{t, x)UU x + g(t, x)UU 1 + b{t, x)U x + e(t, x)U 
+ s(t,x)U 2 + m(t,x)Ui + f(t,x)]dt + UidW{t), < t < T, x£l, 



(69) 



dS {t) = r(t)S (t)dt, < t < T, S (Q) = 1, and 
dS{t) = S(t) [n(t)dt + a(t)dW(t)] , < t < T, S{0) = s>0, 



(70) 
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stock by offering tax alleviation analogous to the size of the investment 7r(-) and to the undertaken 
risk captured by the risky asset's volatility cr(-); in particular, the tax rate is given by 



L(t) = -n(t)a 2 (t)Y 2 (t) + r(t)Y(t). 

Therefore, in accordance with the market dynamics of (|70p . the wealth process Y(-) = Y n,c (-), 
corresponding to the portfolio-consumption pair (ir, C), is the solution of the following SDE 

dY(t) = (1 - it{t))r(t)Y{t)dt + ir{t)Y(t) \fi{t)dt + a(t)dW(t) 
- C(t)dt - L(t)dt 

n(t)(n(t) - r(t))Y(t) + ir(t)a 2 (t)Y 2 (t) - C(t)^dt 

+ Tr(t)a{t)Y(t)dW(t). (71) 

Of course, all the preceding decisions taken by either the agent or the government should depend 
on the information available up to t and not anticipate the future, thus the portfolio strategy 
7T : [0, T] x $7 — > R and the consumption strategy C : [0, T] x Q, — > [0, oo) are assumed to be 
F-progressively measurable processes, and in addition verify the technical integrability condition 
Jq T (C(t) + 7r 2 (t)^dt < oo, almost surely. 

Let us now broach in this market the issue of pricing a contingent claim which depends on 
the volatility of the stock; that is, an J 7 (T)-measurable random variable ^ > 0, satisfying proper 
integrability conditions, of the form £ = p(X(T)), where X(-) is the forward process of (|33p and 
p:lxfl^ [0, oo) is an J-^-measurable random field. This random amount £ represents a liability 
for its seller that has to be covered with the smallest amount Y(0) of initial funds at time t = 
and the right trading strategy ir(-) during the interval [0,T], so that the corresponding wealth 
process Y(t) is positive for all t € [0, T] and at the end of the time-horizon Y(T) = p(X(T)) holds 
without risk. In other words, to hedge and determine the fair price 1^(0) of the contingent claim 
at time t = it suffices to find the solution pair (Y, it) of the BSDE (|7"T|) subject to the previous 
terminal condition such that Y(-) > 0. 

Setting 

Z(t) = Tr(t)a(t)Y(t), < t < T, (72) 
and eliminating the portfolio ir(-) in the drift term of (I71|) . we get the BSDE 

dY(t) = \a(t)Z(t)Y(t) + m(t)Z(t) - C(i)l dt + Z(t)dW(t), < t < T, 
Y(T)=p{X(T)), 

where 

/,(•)- r(-) 



is the well known relative market of risk process. Due to Proposition 15. 1[ its solution is given by 
the pair (Y, Z) of ([32]), where (U,^ u ) is of class C F ([0, T]; h 2 (Q; C 3 (R))) x L|(0,T;C 2 (R)) and 
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solves the BSPDE 
dU 



1 



a 2 (t)U xx + a 2 (t)UU x + a{t)UV u + m(t)a(t)U x - a{t)^ u x 



+ m(t)V u -C(t) 
U(T,x)=p(x), ieR. 



dt + -$ u dW(t), 0<t<T, xG 



(73) 



In turn, Theorem 15.61 provides the solution of this BSPDE via the generalized Cole-Hopf transfor- 
mation (|55D , in terms of a triplet (V, ^ v , ^ V ) G C F ( [0, T] ;L 2 (S7; C 4 (M))) x C F ( [0, T] ;L 2 (S7; C 3 (R))) 
that solves the linear BSPDE 



xL|(0,T;C7 



dV 



■ o\t)V xx + a(t)m(t)F x - <7(t)*£ + m(f)* v + C(i)z;V }> eft 



F(T,x) = g(x), x£R 
for c(-) = 0, subject to V( • , •) > 0, 



v 



-a\t)^ x +[a(t)m(t) + ^(t) 



v 



+ 



a(t) 



m(t) 



a(t) 



+ C(t)x 



o*(t) 



* y - <r(t)*r 



diy(t) on [0,T) 



and (|56p . Then, the hedging portfolio process 7r(-) follows immediately from (|72p . 

Finally, we shall illustrate simple and explicit solutions of the above system in the case of 
constant model coefficients r, fi, a, so as m, and zero consumption, i.e., C(-) = 0. In light of 
Examples 15.101 and 15.111 we postulate strictly positive solutions of the form 



V(t,x) = e aW ^ or V{t,x)=l3e 



W(t) 



for given constants a, (3 > 0, and take a < so that (|55p establishes positive solutions for the 
BSPDE (I73p . respectively. Then, in order for the corresponding 2x2 systems in these examples 
to hold, it suffices to select the coefficients of the model such that 



m(-) = — a 
v ; 2 



and from (|56p 
respectively. 



a 



> 



or 



or 



m(-) = -, 
w 2 



a 



Pi- 



er 



>0, 
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